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ABSTRACT 

Me  discuss  the  nonlinear  eigenvalue  problem 

u"  +  Xf(x,u)  -0  -1  <  x  <  1,  A  £  0 

(P) 

u’ (-1)  -  u*(+1)  -  0,  0  i  u( x )  $  1 

where 

f(x,u)  -  u{ 1  -  u) [u  -  a(x)] 

and 

a(x)  -  <-  [1  -  er(x)  +  hi  e  £  0,  he* 
with  r(-x)  »  -r(x)  and  r'  £  0. 

Por  £  ■  h  ■  0  the  solution  of  Problem  P  is  well  known,  and  every  solution,  except 
u  *  0  and  u  ■  1  is  unstable  with  respect  to  the  corresponding  parabolic  problem. 

we  show  how  the  branch  of  increasing  solutions  changes  as  e  becomes  positive,  and 
acquires  a  bifurcation  point  (X,u)  beyond  which  this  branch  becomes  stable.  If  h 
becomes  nonzero  as  well,  this  bifurcation  point  '<.»  shown  to  break  up. 

As  an  illustration  we  consider  an  example  in  which  the  branch  of  increasing  solutions 
can  be  computed .  Here 


f ( x,  u)  - 


-u  0  $  x  <  a(x) 

1  -  u  a(x)  <  u  $  1 


where  a(x)  is  given  above. 


e  >  0,  h  =  0 


Accession  Por 

NTIS  GRAJfcl 

DTIC  TAB 
Unannounced 

* 

□ 

Justification _ 

By 

Distribution/ 
Availability  Codes 
Avail  and/or 


AMS  (MOS)  Subject  Classifications:  34B15,  35K55  Avail  and 
Key  Words:  nonlinear,  eigenvalue,  bifurcation,  diffusion  equation  Dist  i  Special 
Work  unit  Number  1  (Applied  Analysis)  i  i " 


ft£elft  University  of  Technology,  Delft,  The  Netherlands 
University  of  Leiden,  Leiden,  The  Netherlands. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 


SIGNIFICANCE  AND  EXPLANATION 

Consider  a  population,  consisting  of  three  genotypes t  AA,  Aa  and  aa, 
which  tend  to  disperse  through  a  habitat  ft  and  possess  different  survival 
fitnesses.  Let  u(x,t)  denote  the  fraction  of  allele  A  amongst  the  total 
number  of  alleles  in  the  population  at  a  point  x  in  the  habitat,  at  time 
t.  Then,  according  to  a  model  due  to  Fisher  and  Haldane,  the  evolution  of 
u  is  described  by  the  nonlinear  diffusion  equation 

“  DAu  +  f(x,u)  x  e  ft,  t  >  0 

in  which  D  is  a  positive  constant,  measuring  the  rate  of  dispersal  and  f  a 
function  determined  by  the  relative  fitness  of  the  genotypes.  If  the  fitness 
of  the  heterozygote  Aa  is  lower  than  that  of  either  of  the  homozygotes, 
then  f  is  of  the  form  given  in  the  Abstract.  When  ft  -  (-1,1)  and  we 
assume  no-flux  boundary  conditions  we  arrive  at  Problem  P,  with  X  “  1/D. 

In  this  paper  we  shall  assume  fx  £  0,  and  study  the  set  of 
nondecreasing  equilibrium  solutions,  and  their  stability  properties.  In 
particular,  it  is  shown  that  if  an  equilibrium  solution  is  decreasing  on  some 
Interval,  then  it  must  necessarily  be  unstable. 

Problem  P  contains  three  parameters,  Xt  measuring  the  relative 
Importance  of  dispersal  and  selection,  es  measuring  the  gradients  in  the 
selection  pressure  fx  and  hi  measuring  the  departure  from  the  symmetry 
property 

(*)  f(-x,1  -  u)  •  -f(x,u)  x  6  (-1,1),  u  C  [0,1]  . 

We  find  that  if  (*)  is  satisfied  (h  «  0),  and  fx  %  0,  there  exists 
a  branch  of  increasing  solutions  $(x, X)  such  that 

♦(-x,X)  -  1  -  ♦(x,X)  x  e  [-1,1]  , 

which  is  unstable  for  X  small  and  stable  for  X  large,  the  exchange  of 
stability  being  achieved  at  a  bifurcation  point  on  this  branch.  This 
singularity  disappears  when  h  becomes  nonzero,  and  the  symmetry  relation  is 
no  longer  valid.  Biologically,  this  means  that  for  small  dispersal  rates, 
stable  equilibrium  solutions  exist,  which  mirror  the  selective  advantages  in 
the  habitat. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  NEC,  and  not  with  the  authors  of  this  report. 


ON  A  NONLINEAR  EIGENVALUE  PROBLEM 


OCCURRING  IN  POPULATION  GENETICS 

* 

Ph.  Clement  and  L.A.  Peletier 


1.  INTRODUCTION 


In  a  previous  paper  [8] 
nonlinear  eigenvalue  problem 


the  second  author  considered  the  following 


(P) 


r  u,!  +  Xf(x,u)  “0  'I  <  X  <  1  ,  1  5  0  . 
I  u'(-1)  *  u'(1)  -  0 


where 

(1.1)  f(x,u)  =  u(1-u)(u-a(x)) 

and  a  satisfies 

At  a  €  C1  ([“Ml) 

A2  0  <  a(x)  <  1  for  -1  £  x  S.  1  . 

This  problem  arises  in  population  genetics  as  a  simple  diffusion  model  for  the 
propagation  of  genetic  material  in  a  population,  see  [5,6  j.  The  solutions  of 

are  stationary  solutions  of  the  orohlem 
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(P) 


(D) 


-1  <  x  <  1  ,  t  >  0 
tJO  . 


{u_  ■  u  Xf(x,u) 
t  xx 

ux(t,-1)  -  ux(t,1)  -  0 

Since  u  represents  a  fraction  of  alleles  in  the  population,  we  are  only  inter¬ 
ested  in  solutions  u  taking  values  in  [0,1]  ;  u  ■  0  and  u  ■  1  are  (trivial) 
solutions  of  (P),  and  they  are  both  asymptotically  stable  .  As  a  consequence  of 
the  implicit  function  theorem  one  easily  proves  that  the  branches  { (X,0)  1XGR+} 
and  [  (X,  1)  I XGR+}  do  not  have  a  bifurcation  point  in  R+  x  C([-1,1))  where 
C([-l ,11)  denotes  the  space  of  continuous  functions  on  [-1,1]  equipped  with 
the  supremum  norm. 

When  a  * (x)  **  0  ,  x  €  [-1,1]  ,  the  trivial  solutions  are  the  only  stable 
ones  [1]. 

When  a*(x)  <  0  ,  x  G  [-1,1]  ,  it  was  proved  in  [7, Theorem  4]  and  [5]  that  for 

sufficiently  large  values  of  X  ,  there  does  exist  at  least:  one  non-trivial 

solution  which  is  stable.  (We  prove  in  Appendix  B  that  stable  solutions  have 
to  be  increasing) . 

In  this  paper,  we  are  interested  in  the  question  as  to  how  these  stable 
solutions  emerge  when  a(x)  changes  from  a  constant  to  a  non-constant  function. 

To  be  specific  we  shall  write 

(1.2)  a(x)  -  ]^1-er(x)+h^  ,  e,h  G  B 

where  r(-x)  -  -r(x)  ,  x  G  [-1,1]  . 

% 

In  the  example  given  in  [8],  wjiere  h  »  0  ,  r(x)  -  1  on  (0,1]  and  e  G  [0,1)  , 

it  was  shown  that  the  stable  solutions  are  part  of  continuum  in  B+*C([-1,1]) 

which  extends  all  the  way  back  to  X  »  0  ,  and  that  (X,u)  converges  to  (0,1) 
as  X  4  0  .  In  addition  it  was  shown  in  [8],  that,  under  the  assumptions  Al,  A2, 
for  X  small,  there  exists  a  unique  curve  C  of  nontrivial  unstable  solutions 

*  For  this  choice  of  r  the  function  a  does  not  satisfy  the  condition  Al ,  but 
this  has  no  effect  on  the  result. 
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u(X)  ,  bifurcating  from  (0,a)  where 
+  1 

(1.3)  a  -  |  |  a(x)dx  . 

-1 

Thus,  in  the  example  considered  above,  in  which  f  ha$  the  symmetry 
property: 

(1.4)  f(-x,1-u)  ■  -f(x,u)  x  €  [-1,1]  ,  u  €  [0,1]  , 

(that  is  f  is  antisymmetric  with  respect  to  the  point  (0,J)  in  the 
(x,u)-plane) ,  the  branch  of  nontrivial  solutions  C  extends  to  infinity  and 
eventually  the  solutions  of  this  branch  become  stable.  As  a  consequence  it 
was  shown  that  C  contains  a  bifurcation  point  (X,u)  . 

In  this  paper  we  shall  extend  these  results  to  more  general  functions 
r  ,  and  discuss  the  dependence  of  (X,u)  on  e  . 

If  h  becomes  non  zero  and  f  no  longer  has  the  symmetry  property  (1.4), 
we  prove  some  partial  results  which  suggest  that  the  continuum  C  breaks  up  at 
the  point  (X,u)  .  We  conjecture  on  the  basis  of  these  results  that  when  f 
fails  to  have  the  symmetry  property  (1.4),  then  C  does  not  contain  any  stable 
solution. 

We  conclude  this  paper  with  an  example  of  a  different  type  of  function  f 

r  -u  -0  5  u  <  a(x) 
f(x,u)  -  < 

'  1-u  a(x)  <  u  S  1  , 

where  n(x)  is  given  by  (1.2)  and  r(x)  =  x  .  In  this  example  the  branch  C 
can  be  computed  explicitely,  both  when  li  is  zero  and  when  h  is  not  zero. 

We  gratefully  acknowledge  S'.  Hastings  for  suggesting  this  example. 
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2.  THE  CONTINUUM  C 


Define  the  function  J(u)  * 

Then 


+  1 

/  f (x,u)dx  ,  where  f  is  given  by  (1.1). 
-1 


(2.1)  J(u)  -  O^uS  {0,a,1} 

where  a  is. defined  in  (1.3).  Recall  that  by  A2,  0  <  a  <  1  .In  this  section  we 
shall  show  that  from  the  point  (0,a)  emanates  a  branch  of  nontrivial  solutions  of 
(P).  We  shall  denote  it  by  C.  By  a  solution  of  Problem  (P)  we  shall  mean  a  pair 

(X,u)  G  R  x  C2([-1,1}),  which  satisfies  (P),  where  C2([-1,1])  is  the  space  of  twice 
continuously  differentiable  functions  on  [-1,1]  equipped  with  the  usual  norm. 

THEOREM  1 .  Suppose  a  satisfies  A1  and  A2.  Then  there  exists  a  maximal 
connected  set  of  nontrivial  solutions  C  in  R+  G  C2([-1 , 1 ])  ,  whose  projection 
on  R  is  unbounded ,  with  the  properties: 

(i)  lim(X,u)  «  (0,o)  in  R  *C2([-1,lJ)  . 

(X.uKC 

X+0 

(ii)  There  exists  X^  >  0  such  that  for  X  €  (0,Aq)  there  exists  one  and 
only  one  nontrivial  solution  (X,u)  ,  and  (X,u)  belongs  to  C  . 

(iii)  If  (X,u)  €  C  ,  then  0  <  u(x)  <  1  ,  x  G  [-1,1] 

(iv)  If  a  satisfies  in  addition 

A3  a'(x)  <0  x  G  [-1,1]  , 

a  not  a  constant ,  then  (X*u)  f  C  implies  u'(x)  >  0  ,  x  €  (-1,1)  . 

(u  is  a  "aline” ). 

REMARK.  For  X  6  (0,1^  ,  the  function  X  •*  u(X)  G  C2([-1,1])  is  analytic. 
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PROOF.  To  prove  (i)  and  (ii)  we  shall  rewrite  Problem  (P)  such  that  we  can 


apply  Theorem  A  from  the  Appendix.  Observe  that  if  (X,u)  is  a  solution 

+  1 

of  Problem  (P)  and  X  >  0  ,  then  /  f(x,u(x))dx  =>  0.  Therefore  solutions 

-1 

of  Problem  (P)  with  X  >  0  are  also  solutions  of  Problem  (PI)  and  conversely: 

+  J 

r  -u"  *  x{f(x,u)-i  /  f(x,u(x))dx} 

-1 


(PI) 


\  u'(±1)  -  0 

I  +1 

>  /  f(x,u(x))dx  =  0  . 
-1 


+1 

Let  k  be  such  that  if  g  f.  C(  1-1,1])  with  J  g(x)dx  *  0  and 
+1  -1  +1 
u(x)  :«  / k(x,y)g(y)dy  ,  then  u  satisfies  (i)  j  u(x)dx  =  0  and 

-1  -1  +1 
(ii)  -u"  *  g  ,  x  f.  (-1  ,+1)  and  u'(±1)  **  0  .  Let  3  denote  J  J  u(x)dx 

-1 

and  y  *  3  ~  a  ,  where  a  is  defined  in  (1.3),  and  v  ■  u  -  3  ,  Then  we  can 
rewrite  Problem  (PI)  as 


(P2) 


+  1 

v (x)  **  X  /  k(x,y)f (y,v(y)+o+y)dy 

: 

.  J  f (x,v(x)+a+y)dx  ■  0  . 

-1 


Now  we  define  the  Banach  space 

+  1 

E  |v€C([-1,1])|  Jv(x)dx=oJ  *  R 
-1 

equipped  with  the  norm  l|(v,y)4l  =  max  lv(x)|  +  |y|  ,  and  the  map 

x£[-1 , 1 ] 

F  :=  R  x  e  -*■  E  by 


+  1 

C  v-X  /  k(. ,y)f (y,v(y)+u+y)dy 

P(X;(v,y))  =  <  +1  -1 

^  /  f (x,v (x)t-a+y)dx  . 

-1 
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Plainly  Problem  (P2)  can  be  rewritten  as 

(P3)  F(X;(v,y))  -  0  ,  X  f.  R  ,  (v,y)  e  E  . 

We  are  now  in  a  position  .to  apply  Theorem  A  of  the  Appendix  and  we  shall  prove 
that  the  conditions  a)  to  f)  are  all  satisfied,  a)  follows  from  the  definition 
of  F  and  a'  .  b),  c)  and  f)  with  F  analytic  are  standard.  To  prove  d), 

one  computes  the  partial  Frechet  derivative  of  F  with  respect  to  (v,y) 
at  (0,(0,0))  . 

F(v,y)<0;(0,0))(°,9)  "  {  /+1  f  (x,a)[$(x)+Y]dx  , 

-1  u 

+1 

where  ($,?)  €  E  .  Since  /  f  (x,a)dx  ■  2a(1-a)  0  it  follows  that 

-1  u 

FjVtY)<0,<0,0))  €  Isom(E,E)  .  It  remains  to  prove  e).  First  observe  that  as  a 
consequence  of  (2.1),  F(0,(v,y))  *  0  if  and  only  if  v  -  0  and 
Y  €  {-a,0,1-a}  .  Then  one  observes  that  the  component  of  solutions  in  Rx  E 
which  contains  (0,-a)  (resp.  (0,1-a))  is  exactly  {(X,-a))XGR)  (resp. 

( (X, l-oi)  I XGR}  since  there  is  no  bifurcation  on  these  lines,  thus  the  component 
of  solutions  in  It  x  e  which  contains  (0,0)  does  not  meet  {0}  x  e  ,  and 
condition  e)  is  satisfied.  It  follows  from  Theorem  A  that  Problem  (P3)  possess 
a  maximal  connected  set  of  solutions  in  R+  x  E  which  is  unbounded  in  R+  x  E 
and  such  that 

lira  (X, (v,y))  "  (0,(0,0))  . 

(X, (v,y))GP 
X+0 

o 

It  is  easy  to  verify  that  v  €  C  ([-1,1])  .  Next  we  define 
+  2 

C  :■  {(X,u)GRxC  ([-1  *  1  ])  lu»v+Y+a  with  (X,(v,y))GP)  .  C  is  an  unbounded 

.f. 

maximal  connected  set  in  R  xC([-1,1])  ;  moreover  (X,u)  f.  C  are  nontrivial 
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solutions  of  (P).  Thus  since  the  R'xC([-1,l])  and  the  R*  C2([-1,1]) 

topology  are  the  same  on  the  set  of  solutions  of  (P),  C  is  a  maximal  connected 

set  in  R+  x  C2([-1,1])  ,  such  that  lim(X,u)  =  (0,a)  in  Rx  C2([-1,1])  . 

(A,u)f.C 

X+0 

This  proves  part  (i)  of  Theorem  1 . 


Part  (i  1)  follows  from  f)  of  Theorem  A  .  Next  we  prove  that  (A,u)  G  C  implies 

0  <  u(x)  <  1  ,  x  G  [—1,13  •  It  is  sufficient  to  prove  max  |u(x)~il  <  { 

x€[-1 , 1  ] 

for  (A»u)  G  C  .  First  observe  that  for  X  small,  max  lu(x)-J|  <  {  , 

x€[-1 , 1 ] 

with  (X,u)  €  C  .  Assume  that  for  some  (X,u)  G  C  ,  max  |u(x)~4l  =  J  . 

xGt-1,1] 

Prom  the  boundary  conditions,  it  follows  that  there  exists  £  G  [-1,1]  such 
that  u'(£)  =  0  ,  u(£)  =0  or  1  .  But  the  function  u(x)  b  0  or  a  1  would 
be  a  solution  of  the  initial  value  problem 


•  ~v"  *  Af (x,v) 
(IP)  v(§)  -  u(5) 
-v'(5)  =  0 


as  well  as  the  function  u  which  is  not  constant.  This  contradicts  the  unique¬ 
ness  of  the  initial  value  problem  (IP).  We  have  then  that  (X,u)  G  C  implies 

max  lu(x)*il  -  4  .  One  concludes  by  observing  that  the  function 

xG  [-1,1] 

(A  ,u)  G  C  -*•  max  |u(x)~4l  is  continuous  from  Ex  C([-1,1])  into  R 
x€M,1] 

and  that  C  is  connected  in  ,R  XC([-1,1])  .  This  proves  part  (iii). 

As  a  consequence  we  prove  next  that  Projj,+  C  is  unbounded.  We  have 

0  v(x)  +  Y  ♦  a  <  1  for  all  (X,(v,y))  £  V  .  Since  J+ ^  v(x)dx  =  0  ,  we 

-1 

also  have  0  <  y  +  a  <  1  ,  for  all  (X,(v.y))£P  and  thus  sup  II(v,y)II  < 

(X, (v,y))GP 

hence  ProJ^+  ^  ’-s  founded,  since  V  is  unbounded  in  it+  x  E  . 


Finally  we  prove  (iv).  For  (X,u)  c-  C  ,  we  define  iip(X,u)  to  be  the 
principal  eigenvalue  of  the  problem 

j-  -h"  -  Xfu(.,u)h  ■  yh  -1  <  x  <  1 
h(-1)  *  h(+1)  -  0  . 

For  (X,u)  €  C  and  X  small  it  is  easy  to  see  that  yp(X,u)  >  0  .  But 
PD(X,u)  i  0  when  (X,u)  €  C  .  Indeed  w  ■  u'  satisfies: 


(2.2) 


r  -w"  -  Xf  w  -  Xf  ■  Xu(1-u)(-a')  2  0 
I  u  x 

^  w(-1)  =  w(1)  »  0 


+  1 

Suppose  y  (X,u)  3  0  .  Then  /  f  hdx  should  be  zero  where  h  denotes 
U  -1  X 

the  principal  eigenfunction  chosen  positive  and  normalized  by  max  h(x)  **  1 

x€[-1,1] 

Since  a  is  not  a  constant,  this  is  impossible.  As  in  part  (iii)  we  conclude 
that  uD(X,u)  >0  on  C  by  observing  that  nD  is  a  continuous  function  of 
(X,u)  and  C  is  connected.  Since  the  operator  in  the  left-hand  side  of  (2.2) 
is  coercive,  and  the  right-hand  side  is  nonnegative  (not  identically  zero), 
then  u'(x)  «  w(x)  >0  ,  x  €  (-1,1)  .  This  completes  the  proof  of  Theorem  1. 


3.  THE  SYMMETRIC  CASE 

Throughout  this  section  we  set  h  3  0  .  Thus 

a(x)  -  J(1-er(x))  , 

and  we  assume  that  r  satisfies 

HI  r(-x)  =■  -r(x)  x  C  [-1,1] 

112  r'(x)  0  on  [-1,1]  and  r’(0)  >  0  . 


-8- 


Thus  a  -  J  ,  and  by  Theorem  1  the  continuum  C  "emanates"  from  (0*i)  . 
For  convenience  we  define  v  ■  2u  -  1  .  Then  Problem  (P)  is  equivalent  with 


where 


v"  +  Xg(x,v)  =0  -1  <  x  <  1 

* 

;  v'(-1)  =  v'(1)  "  0  , 


g(x,v)  =  i(1-v  )[v+er(x)3  . 


Suppose  4»*  is  a  solution  of  the  problem 


v"  +  Xg(x,v)  =0  0  «  x  <  1 

v (0)  =  0  v’(1)  -  0 


Then  due  to  the  sysmmetry  property  of  r  given  in  HI,  the  function 


<|>(x)  :> 


-$*(-x)  -1  S  x  <  0 


<J>*(x) 


0  S  x  S  1 


is  a  solution  of  Problem  (P1),  and  <j»(x)  =  iO+T(x))  is  a  solution  of  Problem 
(P).  Since  the  function  g  has  the  property 

g(x,v)  >  gv(x,v)v  0  <  x  S  1  ,  0  s  v  S  1 

the  following  facts  .about  Problem  (P")  are  well-known: 
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PROPOSITION  2.  Let  r  satisfy  hypothesis  HI  and  H2.  Then 

(i)  for  each  X  >  0  ,  Problem  (P")  has  a  unique  positive  solution  $*(X)  ; 

(ii)  f*(X)  is  increasing  and  concave ; 

(iii)  the  map  Xw>  $*(X):  R  ■+  C  ([0,1])  is  analytic 

(iv)  0  £  Xj  <  X2  implies  ♦*(XJ)(x)  <  f*(X2)(x)  for  x  €  (0,1) 

(v)  **(X)  ^  0  as  X  -*•  0  in  02([0,t]). 

(vi)  lim  $*(X)(x)  »  1  for  each,  x  £  (0,1]  . 

•X  -* “ 


They  are  an  easy  consequence  of  the  fact  that  the  principal  eigenvalue 
Pj (X)  of  the  linearized  problem 


-h"  -  Xgv(x,^*(x))h  ■■  ph  0  <  x  <  1 
h(0)  *>  h'(1)  =*  0 


is  positive  for  all  X  >  0  . 
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PROPOSITION  3.  Let  r^  and  r2  satisfy  hypothesis  HI,  H2  and  let  <|>*  and 
be  the  solutions  of  Problem  (P"),  corresponding  to  r1  and  r2  .  Suppose 
>  r2  then  4>*(x)  >  <j>*(x)  for  0  <  x  S  1  . 


PROOF.  Write  f^(x,u)  =  u(1-u) (u-a^(x))  where  =  J(l-er^)  i  =  1,i,  .  Then 


♦f"  +  Xf2(.,**)  =  Xf2(.,**)  -  Xf1  (.,**) 


JeX<j>*(1-$*)  (r2~r1 )  S  0  .  (M)  . 
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Hence  is  a  supersolution  of  the  problem 

v"  +  X  f2(.,v)  -  0 
(  v(0)  “  0  v'(1)  -  0  . 

Since  v(x)  a  0  is  a  subsolution  it  follows  that  there  exists  a  solution  w 
such  that  0  S  v  ?  »  and  because  neither  0  nor  are  solutions, 

0  <  w(x)  <  <fr*(x)  for  0  <  x  2  1  .  Because  Problem  (P")  has  only  one  solution 
by  Proposition  2,  w  *  and  hence  <|>*(x)  <  ♦*(*)  0  <  x  £  1  . 

COROLLARY.  Let  r  satisfy  HI  and  H2  and  let  X  €  H+  be  fixed.  Denote  the  solution 
of  Problem  (P")  by  f*(e)  .  Then  e1  2  c2  implies  2  <J>*(e2)  . 

We  now  return  to  Problem  (P).  We  can  deduce  from  Proposition  2  that 

Cg  {(X,4>(X))U>oJ 

is  an  analytic  curve  of  solutions  with  the  symmetry  property 

(3.1)  +(X)(-x>  -  1  -  *(X)(x) 

such  that 

♦  (X)  -*■  1  as  X  1  0  in  C2([-1,1])  , 
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and 

,  1  if  x  G  (0,1) 

$(A)  (x)  "*■  i  as  A  -*■  «®  . 

1  0  if  x  G  (-1,0) 

This  means,  in  view  of  Theorem  1,  that  C  is  the  subcontinuum  of  C  consisting 

s 

of  the  solutions  with  the  symmetry  property  (3.1).  Observe  that  all  increasing 
solutions  satisfying  (3.1)  belong  to  C 

s 

In  the  following  theorem  we  investigate  the  sign  of  the  eigenvalue  v^(A) 
of  the  linearized  problem 

(  -h"  -  Af  (. ,+<A)  h  -  y(A)h  -1  <  x  <  1 
(EV)  {  U 

1  h'(-1)  -  h'(1)  -  0 

where  (A,$(A))  G  C  .If  they  are  all  positive  for  some  (X,  d>  ( A) )  G  C  ,  we 

S  8 

know  that  $(A)  is  asymptotically  stable  [3]  and  if  one  of  them  is  negative 
that  <J>(A)  is  unstable. 

Let  the  eigenvalues  be  ordered  so  that  vi  <  v2  <  v3  <  *  * •  •  Then  due  to 
the  symmetry  of  f^  ,  for  h  ■  1,2,...  .  We  already  know  that 

Uj  (A)  >  0  for  all  A  >  0  .  Hence  a  Vj  *  >  0  for  all  k  &  2  .  Thus  it 

remains  to  find  the  sign  of  . 

THEOREM  4.  Let  (A,$(A))  G  C  .  Then  there  exists  A .  and\0  (0<A  <A-«*>) 
such  that 

(i)  0  <  A  <  A^  implies  v^(A)  <  0 

(ii)  A;,  <  A  <  »  implies  Vj(A)  >  0  , 

where  v^A)  denotes  the  principal  eigenvalue  of  the  Problem  (EV). 

PROOF,  (i)  This  part  was  proved  in  [ 6 3  - 
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(ii)*We  begin  with  the  observation  that  if  i  ■  1,2  satisfy  HI,  H2  and 

>  r(2)  on  (o,1]  ,  then  the  principal  eigenvalues  vj^(X)  and  vj2^  (X) 
satisfy  the  inequality: 

(3.2)  v{1J(X)  >  v[2)(X)  for  all  X  >  0  . 

This  follows  from  the  fact  that  if  >  r^'  then  by  Proposition  3, 

♦  ^(x)  >  <^2^(x)  for  all  x  €  (0,1)  ,  and  hence  as  follows  from  elementary 
computation  f^1  ^  (x*,$^  ^  (x))  <  f^2^  (x,$^2\x))  for  all  x  £[“1,0)  u  (0,1]  . 

The  variational  characterization  of  the  principal  eigenvalue  of  Problem  EV 
now  gives  the  inequality  (3.2).  In  view  of  the  assumptions  on  r  ,  there  exists 
a  constant  6  >  0  such  that 

xr(x)  £  fix2  for  -1  £  x  <  1  . 

It  follows  from  the  above  observation  that  it  is  sufficient  to  prove  part  (ii) 
for  the  function  r(x)  ■  x  , -1  £  x  S  1  ,  where  the  constant  6  has  been 
absorbed  in  the  factor  e  . 

Let  hj  be  the  principal  eigenfunction  of  Problem  EV  normalized  so 

that  h.  >  0  on  [-1,1]  and  max  h  »  1  .  Then 

1  [-1,1]  1 

(3.3)  v1h1  -  -hj  -  Xfu(.,*(X))b  . 

If  we  multiply  this  equation  by  $'(X)  and  integrate  over  (-1,1)  we  obtain, 
using  the  boundary  conditions  for  h^  and  $  ,  and  dividing  by  X  : 
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+  J  fx(/»*(X)]  hidx“  -2f^1,*(X)(1)J  11,(1)  + 

-1 

+  1 

+  c  |  *a)(i-*a))h1dx  . 

0 

Here  we  have  also  used  the  symmetry  properties  of  f,$  and  h  .  Integrating 
(3.3)  over  (-1,1)  and  dividing  by  X  we  obtain 

+  1 

X-1  v, (X)  |  h, dx  -  -2  fu(.,*(X))  h,  dx  - 

-1  0 
1  1 

-  -2  |  «KX)(1-<KX))h1dx  +  2  J  (2*(X)-1)(<KX)-a)h1dx  . 

0  0 

If  we  multiply  this  equality  by  |  and  add  it  to  (3.4)  we  obtain 

+1 

X"1  v, (X)  J  (♦,+ie)h,dx  -  ~2f(l ,*(X) (1 ))h, (1 )  + 

-1 

1 

+  e  |  ^(X)-1^(X)-a)h,dx  . 

0 

We  assert  that  h, (x)  i  h,(1)  for  0  S  x  S  1  .  Accepting  this  for  the 
moment  we  find  that 


+1 

X"1  v,  (X)  |  (♦'♦ifc)h1dx  £  h,(1)  J-2f(l, <KX)(1))  + 
-1 

1 

+c  |  (2KX)-l)(*(X)-a)dxj  . 


By  Proposition  2,  ♦,(X)(x)  1  as  X 

with  0  <  £  <  1  .  Hence 


uniformly  on  any  interval  [£,1] 


lim  |-2f^l  ,4»(X)(1)y+e 


1 


0 


1 

-  e  (l-a)dx  “  2  +  e  }  xdx)  =  |^,+.  2^  >  ®  ’ 

0  0 


Because  the  coefficient  of  v^(X)  is  positive  for  all  X  >  0  ,  there  exists 
a  X2  >  0  such  that  (X)  >  0  for  all  X  >  X2  • 


It  remains  to  prove  that  h(x)  2  h(1)  for  0  £  x  £  1  .  One  proof  was 
given  in  [8].  Here  we  give  another  proof.  Observe  that  z  ■  hj  is  a  solution  of 
the  problem: 


f  - z "  -  Xf  z  -  v  2  »  Xf  h, 

I  u  1  xu  1 

z(0)  -  z(1)  -  0 

in  which  fx^.,$(X))  -  Je[l-2$(X)]  <  0  .  Note  that  by  the  definition  of 

the  principal  eigenvalue  of  the  operator  on  the  left  hand  side  with  Neumann 
boundary  conditions  is  zero.  Hence  with  Dirichlet  boundary  conditions  it  is 
positive,  which  implies  that  z  <  0  on  (0,1)’,  from  which  the  result  follows. 
This  completes  the  proof  of  the  Theorem  4. 

Since  v^  (X)  <  0  for  all  X  >  0  if  e  -  0  ,  one  would  expect  that 
^1*^2  *  83  E  "*  ©  •  We  shall  show  that  this  is  indeed  true.  Define 

X*  ■  sup 

Clearly  X*  will  depend  on  c.  Thus,  we  shall  write  X*(c)  . 


H 


Vj  <0  on 


[0,X)|  . 


THEOREM  5.  X*(c)  -*•  «  as  e  •*>  0  . 

PROOF.  Suppose  to  the  contrary  that  there  exist  a  sequence  {e^}  ,  e  -*■  0 

as  n  -*•  «°  and  a  constant  H  such  that  X*(e_)  <  M  ,  for  all  n  2  1  . 

i  n 

Then  there  is-  a  subsequence,  still  denoted  by  (e^)  »  an^  X  >  0  such  that 

A*(e  )  -*•  X  as  c  0  . 
in  n 

Since  depends  continuously  on  (X,e),Vj  (X’j'fe^)  «  0  for  all  n  2  1  , 

and  therefore,  by  (3.4) 

1 

f(l  .♦<Xn>0>)  ht(xn(1))  =  ien  |  ♦(Xn)(1-*(Xn))  h^X^dx  , 

0 

where  we  have  written  X  =  X*(e  )  .  Hence 

n  n 

1 

|f(i.*<*nHi>)vxn)(o| s  k  I  k(Vdlt  s  k  • 

0 

By  using  the  continuity  of  $(X)(1)  and  (A)(1)  we  get  by  taking  the  limit 
as  n  -*■  ~  : 

f(l,4>(A)(1))  ht (A)(1)  “  0  . 

But  4>(X)  S  (0,1)  and  because  a(1)  €  (0,1)  ,  f (1  ,«f*(X)(l))  ?  0  .  Since 
h *  (X) .( 1 )  *  0  ,  we  have  hj (A)(1)  ^  0  .  Therefore  we  have  a  contradiction.  This 
completes  the  proof  of  Theorem  5. 
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REMARKS.  1)  If  we  define  X*  »  inf { A>0| Vj (X)>0)  ,  then  0  <  X*  <  X*  <  *  and 
v^X*)  “  v^X*)  "  0  '  s*nce  vi  *s  analytic  i°  *  there  are  at  most  finitely 
many  zeros  of  v^  in  the  interval  ,  and  there  is  at  least  one  zero 

where  changes  sign.  We  denote  by  X  the  first  zero  of  where  changes 

sign.  It  follows  as  in  [7],  that  (X,<j>(X))  is  a  bifurcation  point  from  which 
"emanates"  a.  continuum  of  "nonsymraetric"  solutions. 

2)  We  are  unable  to  prove  that  -^—(X)  <  0  which  would  imply  the  transervality 

condition  needed  for  the  theorem  of  bifurcation  from  simple  eigenvalue.  Since 

dv11 

v  is  analytic  we  only  know  that  there  exists  an  m  £  II  such  that  — — (X)  =  0 

.  2m+1  dX 

for  n  S  2m  and  - -Vv  (X)  >  0  . 

dx2m+1 

3) Concerning  the  nonsymraetric  case,  we  just  mention  the  following  result  which 
suggest  that  the  branch  of  stable  solutions  becomes  disconnected  from 

C  when  h  ^  0  . 


PROPOSITION  6.  There  is  no  6  >  0  such  that  there  exists  a  C1  function 
w:  [-6,6)  x  (X-6.X+6)  -*  C2([-1 , 1 ])  which  satisfies 

i)  w(0,X)  *  C(X)  X  £  (X-6.X+6) 

ii)  w(h,X)  is  a  solution  of  Problem  (P)  for  each  h  £  (-6,6)  . 


PROOF.  If  such  a  6  would  exist,  then  z  ■  -j^(0,X)  would  satisfy  the 
equation: 


-z"  -  Xf  (.,<KX),h=0)z  «=  X4»(X)  (1-4>(X)) 


r  - z  -  AT. 

{ 

1  z’(-1)  -  z’ 


(0=0. 


This  is  obtained  by  differentiating  (P)  with  respect  to  h  ,  at  h  *  0  . 

But  at  (U(X))  ,  the  principal  eigenvalue  v^(X)  vanishes  and  the  corresponding 
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eigenfunction  h,  does  not  change  sign.  Thus  /  X$(1-$)h.dx  +  0  which 

1  -1  +1  _  ' 
contradicts  the  necessary  orthogonality  condition  /  HO-^h.dx  *  0  in  order 

-1 

to  have  existence  of  a  solution  z  .  This  completes  the  proof  of  the  proposition. 
Note  that 


(3.5) 


>,h-o) 


(x)  >  0 


x  <E  [-1,1]  , 


X  <  X  <  X  (next  zero  of  v^(.)  ,  possibly  +  «#)  and  there  is  6  >  0  such  that 


(3.6) 


^X.h-o) 


(x)  <  0 


x  e  [-i,i]  ,  x  -  «  <  x  <  x  . 


The  existence  of  -|^(X,.)  €  C^t-1,1]  follows  from  the  implicit  function 


theorem  and  z 


Tt-  satisfies  (E)  for  X  €  (X-S,,X)  U  (X,X)  for  some 
on  i 


6^  >  0  .  Since  the  right-hand  side  of  (E)  is  positive,  and  the  operator 
on  the  left-hand  side  is  coercive  (v^O)  for  X  €  (X,X)  ,  (3.5)  follows. 
The  inequality  (3.6)  is  a  consequence  of  the  antimaximum  principle  [2]. 


4.  AN  EXAMPLE  IN  WHICH  f(.,u)  IS  PIECEWISE  LINEAR 


f (x,u) 


shall 

assume  that  f 

-u 

0  S  u  <  a(x) 

0 

u  ■  a(x) 

1  — u 

% 

a(x)  <  u  5  1 

where 


(4.1) 


a(x)  :«  ]0“t:x-»h)  0  <  e  <  1  ,  |h|  <  min{e,1-e). 
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Thus  £  is  not  continuous  in  u  but  it  still  has  the  properties 

f(x,0)  -  f(x,1)  -  0 


f(x,u)  <  0  for 

x  e  (- 

■1,1)  ,  0  <  u  <  a(x) 

f(x,u)  >  0  for 

x  €  (■ 

■1,1)  ,  a(x)  <  u  <  1 

We 

consider  the  problem: 

Find  (A, 

,u)  €  R+  x  W2*"(-1,1) 

such 

that 

(P) 

J  -u"  -  \f(.,u) 

,  u'  >  0  on  (-1,1) 

\  u'(-1)  ■  u'(+1)  ■  0 

• 

Thus,  u 

satisfies 

(4.2a) 

-u"  •  -Au  if 

u(x) 

<  a(x) 

(4.2b) 

-u"  »  A(l-u)  if 

u(x) 

>  a(x)  . 

Because  u  is  increasing  and  a  is  decreasing  there  exist  a  unique 
£  €  (-1,1)  such  that  u(x)  <  a(x)  on  (-1,£)  and  u(x)  >  a(x)  on  (£,1] 
Thus  by  (4.2a,b)  and  the  boundary  conditions 

u(x)  -a  cosh  \A(x+1)  if-  -1  S  x  S  { 

u(x)  ®  1  -  B  cosh  VX(x-l)  if  £S  x  S  1  , 

where  a, 3  G  R  . 
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I.  THE  SYMMETRIC  CASE  (h=0) 

Equation  (4.6)  now  reduces  to. 

(4.7)  ?|<5>(--\Af;)-r.4!(2\/A)|  =  0  . 
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Plainly,  £  =  0  is  a  solution  for  all  X  >  0  .  The  corresponding  values  of  a 
and  8  can  be  computed  from  (4.3a).  They  are 

o  *  B  **  1/(2  cos  h  V*) 


Thus,  there  exists  a  continuous  branch  of  symmetric  solutions. 
Other  solutions  of  (4.7)  must  satisfy 


(4.8) 


P(X,£) 


»(2VA  O  m  . 
*(2v/x) 


where  we  have  divided  (4.7)  by  £4>(2^)  .  Because  P  has  the  properties 

(i)  P(X,t)  »  P(X,-t)  for  all  t  €  R  ; 

(ii)  P(X,0)  «  1/4>(2vA)  ; 

(iii)  t  |~(\,t)  >  0  for  all  t  €  IN  (0) 

it  follows  that  (4.8)  has  exactly  two  solutions  ±  £(X,e)  if  and  only  if 
c  >  1/*(2\/X)  ,  i.e.  when  X  >  X^(e)  where  *(2^”(e))  ■  1/e  . 


Fig.  1  The  rase  X  >  Xy 
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J^r|3— 


The  corresponding  values  of  a  and  -8  can  again  be  determined  from  (4.3a). 
Thus,  summarizing  we  have  found 

1.  For  each  X  >  0  ,  there  exists  exactly  one  symmetric  solution,  characterized 
by  5*0. 

2.  For  X  >  Xq(e)  there  exist  exactly  three  solutions:  the  •symmetric  one 
(5*0)  and  two  non-symmetric  ones  characterized  by  5  ■  ±  5(X)  .  It  is  readily 
seen  from  Fig.  1  that 


X  (e)  -*■  { 
0  l 


0  as  e  -*■  1 


a#  as  e  -*•  0  . 


II.  THE  ASYMMETRIC  CASE  (h|<0) 

It  follows  from  (4.5)  that  5  ■  0  implies  X  ■  0  .  Thus  when  X  >  0  ,  we 
may  divide  by  5  ,  and  solutions  must  satisfy  the  equation 


(4.9)  P(X,5)  -  e  -  |  . 


We  consider  the  case  h  >  0  . 

One  verifies  that  for  each  X  >  0  ,  the  function  5  P(5>X)  is  strictly 

convex  and  even.  Thus  P(0,X)  ■  1/$(2vft)  is  the  minimum  of  P(5»X)  . 

The  function  £  -*•  e  -  (h/5)  is  strictly  increasing  for  £  €  [-1,0)  U  (0,1] 
Since  by  (4.1) 


e  +  h  <  P(-1,X)  -  1 


lim  e-  £  >  P(0,X)  -  UH2VX) 

5  to 
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there  exists  for  eacfh  X  >  0  and  h  €  (0,1-e)  exactly  one  solution 


Sj(X,h)  €  [-1,0)  of  equation  (4.9).  Since 

£  ^  ‘  0  for  «  <  0 

the  implicit  function  theorem  implies  that  as  a  function  of  X  is  smooth 

and  strictly  decreasing. 

For  £  €  (0,1]  ,  both  the  left-hand  side  and  the  right-hand  side  of  (4.9) 
are  increasing  in  5  ,  but  5  ■*  P(£,X)  is  strictly  convex  and  £  -*■  c  -  (h/O 
is  strictly  concave.  Moreover,  we  have 

lim  P(£,X)  «  1/*(2VX)  >  lim  (e-  j!) 

£40  £+0  * 


and 


P(1,X)  ■  1  >  e  -  h  by  (4.1). 

Thus  equation  (4.9)  has  zero,  one  or  two  solutions. 

Observe  that  P  as  a  function  of  X  is  strictly  increasing,  because 


P(S,A)  -  P(£,X)  |f(2v^C)-F(2VX)} 


where  F(t)  -  t$’(t)/#(t)  and. 


F*(t)  ■  4>  A  )  “ 


“  J(t«t>)  ^  ^sinh  2t-2t^  >  0  for  .  t  >  0 
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For  X  <  Xq(c)  there  is  no  solution  since 

inf  P(£,\)  >  t  >  c  -  h  ■  sup,  (e-j)  . 

C€(0,1]  ^(0,1]  * 

Because  P(£,X)  ■»  0  as  X  -*>  »  and  5  £  (0,1)  fixed,  there  exist  two  solutions 
of  (4.9)  for  X  large.  Remembering  that  3P(£,.)/3X  <  0  we  conclude  that 
there  exists  a  X*  £  (Xq,«°)  such  that  there  exists  no  solution  of  (4.8) 
if  X  <  X*  ,  exactly  one  of  X  =  X*  and  exactly  two  of  X  >  X*  .  Since 

£[?(«, -  ^P(C.X)  0 

these  solutions  can  be  smoothly  parametrized  by  C  f  (0,1)  . 


For  h  <  0  the  situation  is  entirely  analogous.  This  case  can  be  reduced 


to  the  former  by  substituting  h  -*■  -h  and  £  -*■  -£  in  (4.9). 


APPENDIX  A 


In  this  Appendix  we  state  and  prove  a  version  of  a  global  function  theorem  which 
was  used  in  section  2.  This  variant  is  a  slight  extension  of  Theorem  6.2  of  [9]. 
For  the  sake  of  completeness  we  give  a  oroof  here. 

THEOREM  A.  Let  E  be  a  real  Banach  space  and  let  F:  Ex  E  +  E  satisfy 

a)  F(0,0)  =  0 

b)  F  is  continuous 

c)  G:  RX  E  +  E  defined  by  G(A,w)  =  w  -  F(A,w)  maps  bounded  subsets  of 
Rx  e  into  relatively  compact  subsets  of  E  . 

d)  F  possess  a  partial  Frdchet  derivative  with  respect  to  w  at  (0,0)  : 

dF  (0,0)  and  dF  (0,0)  €  Isom(E,E)  . 
w  w 

e)  If  J  denotes  the  set  of  solutions  of 

(1)  F(A,w)  -  0  in  E*E 

and  if  D  denotes  the  maximal  connected  subset  of  J  which  contains  (0,0)  , 
then 

D  0  ({0}*E)  -  (0,0)  . 

Then  D  is  unbounded  in  R+  *  E  and  in  R  *  E . 

f)  If  moreover  F  €  (^(R’EjE)  ,  then  there  exists  6  >  0  ,  >0  and  a 

function  X  -*■  z(A)  F.  E  such  that  (A,u)  €  J  and  |A|  >  X^  ,  Hull  <  <5  , 

imply  u  *=  z-(X)  .  The  function  z  is  Cm  (resp.  analytic)  whenever  F  is 
Cm  (rasp,  analytic). 


PROOF.  If  the  conclusion  of  the  theorem  does  not  hold,  then  either 

C+  :■  C  fl  (R+*E)  or  C  :■  C  (1  (E  *E)  is  bounded.  Assume  C+  bounded. 

It  follows  from  b)  and  c)  that  C*  is  a  compact  metric  space  under  the  induced 

» 

topology  from  Ex  E  .  It  follows  from  d)  that  (0,0)  is  an  isolated  solution 

of  (1)  in  {0}  x  E  .  Since  C+  is  compact  it  follows  from  e)  that  there 

exists  a  6— neighbourhood  of  C+  in  ]R  *  E  denoted  by  U  such 
U  0  ({0}flE)  fl  J  ■  (0,0)  .  Let  K.:=  U  fl  J  where  U  denotes  the  closure  of  U 
in  Ex  E  .  Since  IJ  is  bounded,  U  is  bounded  and  thus  U  fl  J  is  totally 
bounded.  U  and  J  are  closed,  thus  K  is  compact  metric.  By  construction 
3U  fl  C+  =  0  .  Hence  by  [9],  there  exist  disjoint  compact  subsets 
A,B  e  K  such  that  C+  c  A  ,  3U  fl  J  cz  B  and  K  =  A  U  B  .  Let  0  be  a 
p-neighbourhood  of  A  where  p  is  less  than  the  distance  from  A  to  B 
and  less  than  the  distance  from  A  to  311  .  Note  that  these  distances 
are  positive  since.  A  is  compact  and  B,  3U  are  closed,  A  fl  B  m  0  and 

A  fl  3U  ■  0  .  Hence  C+  c0  ,  30  fl  J  -  0  0  is  bounded  and  0  fl  ({0}*E)  =  (0,0) 

0  fl  ({X}*E)  =  0  for  X  large  enough.  Define  0^  :=  (ufE  I  (X,u)G£)}  .  From 
b)  and  c) ,  it  follows  that  the  Leray-Schauder  degree  of  F  is  well-defined  on 
O^:  d(F,0^)  .  For  X  large  enough,  d(F,0^)  =  0  since  0^  =  0  .  For  X  =  0 

d(F,0rt)  *  index(F(0, . )  ,0)  «*  index(dFv.(0,0)  ,0)  f*  0  since  dF  (0,0)  G  Isom(E,E) 
But  this  contradicts  the  homotopy-invariance  of  the  Leray-Schauder  degree  (see 
Lemma  1.8  of  [9]).  Thus  C+  is  unbounded  in  H+  x  e  .  Similarly  C  -is 
unbounded  in  R  x  E  . 

The  last  assertion  follows  directly  from  the  local  implicit  function  theorem. 
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APPENDIX  B 


The  aim  of  this  Appendix  is  to  prove  that  stable  stationary  solutions  of  (D) 
are  increasing  when  a  is  nonincreasing.  This  generalizes  an  earlier  result 
of  Chafee  [1]  .  It  is  known  that  if  (X,u)  is  a  stable  stationary  solution 
of  (D) ,  then  the  principal  eigenvalue  v  ^  of  •  (EvN) 

t  -h"  f  (.  ,u)h  ■  vh 
(EvN)  <  U 

1  h* (-1)  -  hf(1)  -  0 

is  nonnegative.  Thus  the  principal  eigenvalue  of  (EVD) 


(EvD) 


-h"  -  fu(.,u)h  -  vDh 


f  “n  "  V 

l  h(-1)  -  h( 


is  strictly  positive:  vD  ^  3  v^  2  0  . 
PROPOSITION.  Let  a  eatiafy  At  ,A2  and  let 

A3'  a'(x)  SO  x  C  [-1,1]  . 


Let  (X,u)  be  a  stable  solution  of  (P)  satisfying  0  S  u(x)  SI  ,  x  €  [-1,*1]  » 
X  >  0  .  Then  either  u  is  identically  constant,  u  «  0  ,  u  ■  1  or,  when  a 

is  not  a  constant ,  u  is  strictly  increasing. 


PROOF,  v  *  u'  satisfies 


|  -v"  -  Xfu(.»u)v  ■  Xf^  ■  Xu(1~u)(-a')  i'  0 
^  v<-1)  »  v(1)  «  0  . 


-27 


MV  i  iihmvi  jiiiiyiM  n  >  i 


Since  vn  >  0  %  the  operator  on  the  left-hand  side  is  coercive.  If 
a  constant,  then  it  follows  from  the  maximum  principle  that  0  <  u(x) 
x  £  [-1,1]  and  thus  the  right-hand  side  is  nonnegative  (^0).  Then 
u'(x)  *  v(x)  >0  on  (-1,1)  .  If  u  is  a  constant  it  can  only  be  a 
solution  if  a  is  not  a  constant,  and, if  a  is  a  constant,  and  u  - 
easily  seen  to  be  unstable. 


u  is  not 
<  1  , 

trivial 
a  ,  it  is 
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